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Time Variability of the Neutral Iron Lines from the Sgr B2 Region and 
its Implication of a Past Outburst of Sgr A* 
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Abstract 

We investigate long-term X-ray behaviors from the Sgr B2 complex using archival data of the X-ray 
satellites Suzaku, XMM-Newton, Chandra and ASCA. The observed region of the Sgr B2 complex includes 
two prominent spots in the Fel K-a line at 6.40 keV, a giant molecular cloud M 0.66—0.02 known as the 
"Sgr B2 cloud" and an unusual X-ray source G 0.570— 0.018. Although these 6.40 keV spots have spatial 
extensions of a few pc scale, the morphology and flux of the 6.40 keV line has been time variable for 10 
years, in contrast to the constant flux of the FeXXV-Ka line at 6.67 keV in the Galactic diffuse X-ray 
emission. This time variation is mostly due to M 0.66— 0.02; the 6.40 keV line flux declined in 2001 and 
decreased to 60% in the time span 1994-2005. The other spot G 0.570—0.018 is found to be conspicuous 
only in the Chandra observation in 2000. From the long-term time variability (~10 years) of the Sgr B2 
complex, we infer that the Galactic Center black hole Sgr A* was X-ray bright in the past 300 year and 
exhibited a time variability with a period of a few years. 

Key words: Galaxyrcenter — ISMrclouds — ISM:individual (Sagittarius B) — X-rays:individual 
(Sagittarius B) — X-rays:ISM 



1. Introduction 

Our Galactic Center (GC) black hole, Sgr A*, has a 
mass of 4 x 10 6 M Q (Eisenhauer et al. 2005; Ghez et al. 
2005), but is quiescent with an X-ray luminosity of only 
10 33-34 ergs s _1 (Baganoff et al. 2001), which is many 
orders of magnitude smaller than that of canonical AGNs. 
Sgr A* has been also found to exhibit frequent X-ray flares 
(e.g., Baganoff et al. 2001). This X-ray flare activity is 
not huge, however, with a typical duration of 10 ks and 
an amplitude that is only a factor of ten above the non- 
flare phase. Therefore, even at the flare peak, the X- 
ray luminosity is ~ 10 6 times lower than the Eddington 
luminosity of this supcrmassivc black hole. 

The 6.40 keV line features in molecular clouds suggest 
an active phase in the past. The Sgr B complex of giant 
molecular clouds emit peculiar X-rays with a strong K- 
shell transition line at 6.40 keV from neutral iron (equiv- 
alent width of 1-2 keV) and a deep absorption edge at 
7.1 keV (equivalent 7V H of 10 23 " 24 ) (Koyama et al. 1996). 
The 6.40 keV line flux and the spectrum necessitate the 
existence of strong X-ray sources irradiating the clouds. 
The detailed X-ray morphology of the neutral iron line in 
the Sgr B2 cloud has been investigated by Chandra and it 
was found that the external source is likely to be located 
in the GC direction (Murakami et al. 2001b). However, no 
such source has been found, neither toward the GC, nor 
in regions near the GC. Based on these results, Koyama 
et al. (1996), Murakami et al. (2001b) and Koyama et al. 
(2007b) proposed that the Sgr B2 complex contains X-ray 
reflection nebulae (XRNc) irradiated by Sgr A*, which 
must have been six orders of magnitude X-ray brighter 



about 300 years ago, which is the light travel time be- 
tween Sgr B2 and Sgr A*. Other molecular clouds such 
as M 0.74-0.09, Sgr Bl, Sgr C and M 0.11-0.08 also have 
intense 6.40 keV lines (Koyama et al. 2007b, Nobukawa 
et al. 2008, Murakami et al. 2001a, Yusef-Zadeh et al. 
2007). Some, if not all, of them, arc also likely to be XRNc 
(Koyama et al. 2007b, Nobukawa et al. 2008, Murakami 
et al. 2001a). 

Since the 6.40 keV line can be also produced by high- 
energy electrons bombarding the cloud (Tatischeff 2003), 
Predehl et al. (2003) and Yusef-Zadeh et al. (2007) ar- 
gued against the XRN scenario and proposed electron 
brcmsstrahlung instead. One way to distinguish between 
these two possible scenarios (i.e., irradiation by X-rays or 
electrons (chargcd-particles)) would be to detect the time 
variability of the 6.40 keV line. 

Muno et al. (2007) discovered the time variability of the 
6.40 keV line clumps in the molecular cloud M0. 11— 0.11. 
Koyama et al. (2008a) also found changes in the flux and 
morphology of the 6.40 keV line of the Sgr B2 complex 
between the Chandra and Suzaku observations in 2000 
and 2005, respectively 

Since variability over a few years from a pc-scale ob- 
ject requires light-speed communication, charged particles 
are unlikely to be the origin of the 6.40 keV line. This 
time variability, on the other hand, supports the XRN 
scenario; X-rays from external sources irradiate molecular 
clouds and produce the fluorescence iron line at 6.40 keV 
and Thompson scattering continuum emission in the same 
manner as an X-ray echo. Using the time variable X-ray 
echo, we can infer the light curve of Sgr A* in the past. 

This paper reports a comprehensive X-ray study based 
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on all the available archival data of the Sgr B2 com- 
plex (i.e., the ASCA, Chandra, XMM-Newton and Suzaku 
data). We derived a 10-year-scale light curve for the 
6.40 keV line. We found the time variability of the 
6.40 keV line from the Sgr B2 complex by using the 
time-constant 6.67 keV line from the GC diffuse X-rays 
(GCDX) as a reference. In this paper, the distance to 
the Sgr B2 complex is assumed to be 8.5 kpc (Reid et al. 
1988). 

2. Observations and Data Reduction 

We used four deep exposure and two short survey obser- 
vations as listed in table 1. However, for imaging analysis 
we used only the deep exposure observations. 

2.1. Suzaku 

The Sgr B2 complex was observed with the XIS on 10- 
12 October 2005. The XIS consists of four sets of X-ray 
CCD camera systems (XIS0, 1, 2, and 3) placed on the 
focal planes of four X-ray telescopes (XRTs) on board the 
Suzaku satellite. XIS0, 2 and 3 have front-illuminated 
(FI) CCDs, whereas XIS1 has a back-illuminated (BI) 
CCD. Detailed descriptions of the Suzaku satellite, XRT 
and XIS are found in Mitsuda et al. (2007), Serlemitosos 
et al. (2007) and Koyama et al. (2007a). The XIS observa- 
tion was made in the normal mode. The effective exposure 
time after removing the epoch of the low earth elevation 
angle (ELV <5° ) and the South Atlantic Anomaly was 
about 89 ks. We analyzed the data using the software 
package HEASoft 6.2. The XIS gain was fine-tuned using 
the same procedures described in Koyama et al. (2007b). 
Using point sources in the field, we re-registered the ab- 
solute astrometry by shifting (-0.°0001, -0.°0062) in the 
(I, b) coordinate, referring the Chandra point-source po- 
sitions (Koyama et al. 2007b). 

2.2. XMM-Newton 

The X-ray data were obtained using the European 
Photon Imaging Camera (EPIC) (Struder et al. 2001; 
Turner et al. 2001) onboard the XMM-Newton on 1 April 
2001 and 4 September 2004. The observations were per- 
formed in extended full-frame mode. The data were ana- 
lyzed using Science Analysis Software (SAS 7.1.0). Event 
files for both the PN and the Metal Oxide Semiconductor 
(MOS) detectors were produced using the epproc and 
emproc tasks of SAS, respectively. The event files were 
screened for high particle-background periods. In our 
analysis, we dealt only with events corresponding to pat- 
terns 0-4 for the PN and 0-12 for the MOS instru- 
ments. For the observation in 2004, we made the to- 
tal band images (0.3-10 keV) for each detector and de- 
tected point sources (93 point sources from MOS1, 97 
from MOS2, and 111 from PN). 27 sources were detected 
with all three detectors within positional uncertainties of 
10 arcsec, the sizes of PSF of EPIC at 6 keV, and 19 
sources have possible counterparts with catalogued point 
sources of Chandra (Muno et al. 2003). The mean po- 
sition difference of these point sources relative to the 



Chandra positions is (— 0.°0006, 0.°0003) in celestial coor- 
dinates. We fine-tuned the XMM-Newton frame by shift- 
ing (Act, AS) = (0.°0006, -0.°0003). 

2. 3. Chandra 

Sgr B2 was observed using the Advanced CCD Imaging 
Spectrometer (ACIS-I) onboard the Chandra observatory 
(Weisskopf et al. 2002) on 29 March 2000 and 16 July 
2001. We used the event files provided by standard 
pipeline processing. Only the grade 0, 2, 3, 4, and 6 events 
were used in the analysis. Each CCD chip subtends an 
8.3-arcmin square in the sky, while the pixel size is 0.5 
arcsec. The on-axis spatial resolution is 0.5-arcsec (full 
width at half-maximum (FWHM)). Images, spectra, an- 
cillary files and response matrices were created using the 
software CI AO v3.4. The absolute positional accuracy is 
0.6 arcsec (Weisskopf et al. 2003). 

2.4. ASCA 

The ASCA observations of Sgr B2 were performed on 
22 and 24 September 1994 for ~ 80ks in total. We used 
the event files provided by standard pipeline processing. 
We merged the two observation data since the observa- 
tion dates and aim points are almost the same. Since 
the CIS data did not have sufficient energy resolution 
(~ 470 eV at 5.9 keV) to separate the iron lines, we used 
only SIS data for the analysis. Unlike Suzaku and XMM- 
Newton, no available point sources to use as a reference 
for Chandra for the absolute astrometry correction were 
found. Therefore, the astrometric uncertainty of SIS im- 
age remains about 40 arcsec 1 . We used the software pack- 
age HEASoft 6.2 in the analysis. 

3. Analysis 

3. 1 . Spectra of the Sgr B2 region 

Figure 1 shows relevant diffuse sources in the Sgr 
B2 complex overlaid on the XMM-Newton 6-7 keV im- 
age. Since the FOVs and pointing positions of the 
Suzaku, Chandra, XMM-Newton and ASCA observations 
are slightly different from each other, we selected a region 
in which all the observed areas overlapped (referred to 
hereafter as the Sgr B2 region). The Sgr B2 region (the 
overlapping region) is indicated by by the solid squares in 
figure 1 . The subsequent data analysis and discussion are 
based on the sources in the Sgr B2 region. Hereafter, we 
refer to the Sgr B2 cloud as M 0.66-0.02. G 0.61+0.01 is 
a new source with a strong 6.67 keV line, hence possibly a 
new young SNR (Koyama et al. 2007b). G 0.570-0.0018 is 
an unusual diffuse source discovered by Senda et al. 2002. 
M 0.74—0.09 is also a new 6.40 keV source, but it is not 
discussed in this paper, because it is outside of the Sgr B2 
region. 

We extracted X-ray spectra from the Sgr B2 region and 
subtracted the off-plane blank-sky regions. Thus, the cos- 
mic X-ray background (CXB) and non- X-ray background 
(NXB) were subtracted but the GCDX were not sub- 

1 See http://heasarc.gsfc.nasa.gov/docs/asca/cal_probs.html 
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Table 1. Observation Log 



Observatory 


UDS. id 


Date 


Exposure Time 






lyyyy/ mm/ aaj 


(sec) 


ASCA 


52006000 


1994/09/22 


58,395 


ASCA 


52006001 


1994/09/24 


21,500 


Chandra 


944 


2000/03/29 


98,629 


XMM-Ncwton 


0112971501 


2001/04/01 


12,649 


Chandra 


2280 


2001/07/16 


10,423 


XMM-Ncwton 


0203930101 


2004/09/04 


42,282 


Suzaku 


100037060 


2005/10/10 


89,072 





Fig. 1. XMM-Newton image in the 6-7 keV band. The solid 
line indicates the selected area of the Sgr B2 complex (see 
text). The dashed circles are the individual diffuse sources 
discussed in this paper. M 0.66—0.02 is the Sgr B2 cloud. 
G 0.61+0.01 is a new source with a strong 6.67 keV line, pos- 
sibly a new young SNR. G 0.570— 0.0018 is an unusual diffuse 
source. 

tracted. As the off-plane blank-sky we selected the north 
ecliptic pole data for XIS, distributed blank-sky databases 
for MOS and ACIS, and Lockman hole data for PN and 
SIS. 

The blank-sky-subtracted spectra are shown in figure 
2. These spectra exhibit continuum emissions and pro- 
nounced peaks due to Fel Ka (6.40 keV), FeXXV Ka 
(6.67 keV), and the composite of FeXXVI Ka (6.96 keV) 
and Fel K/3(7.06 keV). The continuum emission is com- 
posed of: 1) the continuum of the GCDX (Koyama et al. 
2007c), 2) the continuum (thermal brcmsstrahlung) from 
G 0.61+0.01, a new SNR in this region (Koyama et al. 
2007b), and 3) a non-thermal power-law continuum with 
a deep iron K-edge related to the 6.40 keV line (Koyama 
et al. 1996). 

It is practically impossible to resolve these three con- 
tinuums; we therefore represented these three continuum 
by a single power-law continuum with absorption (the 
wabs*power-law model in the xspec code). We then added 
four K-shell lines by Gaussian functions. These arc the 
FeXXV Ka (~6.67 keV) and FeXXVI Ka (~6.96 keV) 
lines due to the GCDX, and Fel Ka (6.40 keV) and Fel 
K/3(7.06 keV) lines due mainly to the Sgr B2 region. Thus, 



the model spectrum in the xspec code is given as; 

wabs*power-law + 4 Gaussians (1) 

3.1.1. Suzaku 

The spectra of the Sgr B2 region and background sub- 
tractions of the blank sky for XIS 0,1,2 and 3 were per- 
formed separately, however all the FI-CCD spectra (XIS0, 
1 and 3) were combined, because their response functions 
are almost identical with each other. We fitted the FI- 
CCD and BI-CCD spectra simultaneously in the energy 
band of 5-8 keV with the model given by equation (1). All 
the parameters were free under the following constraints: 
the line widths were fixed to be eV except that of the 
FeXXV Ka line (Koyama et al. 2007c), and the energy in- 
terval between Fel Ka (6400 eV) and K/3 (7058 cV) was 
fixed at the theoretical value of 658 cV (Kaastra & Mewe 
1993). The best-fit spectra for the FI and BI are shown 
in figures 2(a) and 2(b) respectively, while the best-fit pa- 
rameters are listed in table 2. 

Although the K-shell lines from highly ionized iron, the 
FeXXVI Ka and FeXXV Ka lines are due to the largely 
extended GCDX, the best-fit center energy of the FeXXVI 
Ka line (6665 eV) and the flux ratio of FeXXVI Ka to 
FeXXV Ka (0.30) are slightly smaller than those (6680 eV 
and 0.34, respectively) found in the GC region (Koyama 
et al. 2007c). This may be due to the contamination 
of a strong FeXXV Ka line of a new SNR G 0.61+0.01 
(Koyama et al. 2007b). 

When fitting the XMM-Newton, Chandra and ASCA 
spectra, we fixed the line centroids and the Ka flux ratio 
of Fe XXVI and Fe XXV to the best-fit Suzaku values given 
in table 2, because the statistics and energy response of 
these satellites were limited compared to Suzaku. 

3.1.2. XMM-Newton 

We fitted the MOS and PN spectra separately, where 
the line center energies and the Ka line flux ratio of 
FeXXVI and FeXXV were fixed to those of the Suzaku 
best-fit values (see Section 3.1.1). This fit failed in re- 
producing the line and edge structures in the PN spectra, 
indicating a systematic gain shift; we thus added a gain 
offset as a free parameter. The fit was then acceptable as 
Table 3 shows. The best-fit spectra are given in figures 
2(c)-(f). 

3.1.3. Chandra 

We fitted the Chandra spectra using the same method 
as the XMM-Newton case. The spectrum in 2000 was well 



(a) Suzaku XIS-FI 2005 
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(b) Suzaku XIS-BI 2005 
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Fig. 2. The blank-sky-subtracted spectra of the Sgr B2 region, with the Suzaku XIS FI sensors (a) and BI sensor (b) in 2005, 
XMM-Newton MOS (c) and PN (d) in 2004, (e) and (f) in 2001, Chandra ACIS-I (g) in 2001 and (h) in 2000, and ASCA SIS (i) in 
1994. The solid lines indicate the best-fit models (see text). 



fitted by adding a gain offset of ~ —20 eV. The best-fit 
spectra and parameters are given in figure 2(g) and table 
3. In the 2001 spectra, however, the best-fit photon in- 
dex was unrealistic with T ~ —2.2. If we fixed T = 2.6, 
the best-fit Suzaku value, then there were large residuals 
above the 7 keV band. We compared the count rates in 
the 10.0-12.5 keV band between the source and blank-sky 
spectra, and found that the high-energy band count-rate 
of the blank-sky was 2.4 times smaller than that of the 
source spectra. Thus, the strange behaviors are likely due 
to unstable NXB. We therefore subtracted the blank-sky 
spectra by first multiplying them by a factor of 2.4. We 
then obtained a good fit with reasonable parameters. The 
best- fit spectra and parameters are shown in figure 2(h) 
and table 3, respectively. It could be argued that the 
background subtraction described above is artificial, and 
hence may cause large systematic errors in the line flux. 
However, the line flux errors for the different background 
subtractions are within 10%. More details concerning the 
line flux dependence on the continuum shape and/or back- 
ground subtractions are given in Section 3.2.2, for the case 
of the XMM-Newton (table 5). 
3.1.4. ASCA 

The fitting procedures for the ASCA spectra were the 
same as those for XMM-Newton and Chandra spectra. 



The ASCA spectra require a large systematic gain offset 
about 80 eV to adjust the iron line centroids. This large 
offset is caused by chip-to-chip variation of the gain 2 . The 
best-fit spectra and parameters are shown in figure 2(i) 
and in table 3. 

3. 2. Systematic Errors oj the Line Flux 

The estimation of the line flux depends on the under- 
lying continuum shape and the background subtraction. 
We thus examined the possible systematic error cause by 
the assumed continuum and background subtraction. 
3.2.1. Dependence of the Line Flux on the Continuum 
Shape 

We fixed either T (Case 2) or 7V H (Case 3) or both 
(Case 4) to those of the best-fit Suzaku results and fit- 
ted the XMM-Newton, Chandra and ASCA spectra. For 
comparison, we call fitting with free T and iVu as Case 1. 
The best-fit parameters of these fittings (Cases 2, 3, and 
4) are shown in table 4. 

From table 4, we find that the best-fit fluxes of the 
6.40 keV and 6.67 keV lines for the four fitting cases are 
consistent within ±5% with each other, which is smaller 
than the statistical errors. We therefore ignore these sys- 

2 See http://heasarc.gsfc.nasa.gov/docs/asca/4ccd.html 
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Table 3. The XMM-Newton, Chandra and ASCA fitting results for the Sgr B2 region 



Observatory 


XMM-Newton 


XMM-Newton 


Chandra 


ASCA 


Detector 


MOS(2004) 


PN(2004) 


MOS(2001) PN(2001) 


ACIS(2001) 


ACIS(2000) 


SIS(1994) 








— Continuum — 








r 

rpow 


4.8i°J 


1 +0.4 

Z ' i -0.4 

1 fi +0 ' 2 
^•"-O.l 

7 o +0 - 2 

• - u -0.2 


0.5 Q ^ 1.5 Q Q 

0.431 0.93! 
S.StH 9At° 3 3 
— Neutral iron lines — 


i-0±i;f 
7 9+0.5 

' - z -0.5 


1 o +0 - 5 

7+ - 4 
u - ' -0.5 

6 7+ - 2 


0.7+" 

— U. 1 

1 l+°' 7 
5.8±g;| 


^640 
F706 


J - uo -0.08 

n i*?+ 07 

u - lo -0.07 

n qi+ 004 

u - al -0.08 


1 06+ 05 

1 - UD -0.07 
u - 10 -0.15 

u -° y -0.06 


1 co+0.15 -1 97 +0.11 
1 - dz -0.16 1 ' z '-0.12 

U - 1D -0.14 u - 1( -0.09 

— Highly ionized iron lines — 

1 99+O.I5 n 01 +0.11 
4 - zz -0.16 u - Si -o.ii 


1 07+O.26 
i - o '-0.26 

10+ 032 

1 1 4+0-27 
1 - 1 ^-0.27 


1 7-, +0.08 
1 - ' -0.08 

17+0- 09 

u - 1 '-0.04 

1 00+ 08 

1 - uu -0.07 


1 fiS+ 018 

22+ - 15 

u - zz -0.13 

1 20+ 018 

1 - zu -0.18 


X 2 /dof 


0(fix) 
88.9/93 


21+J 
185.6/169 


O(fix) 35+^ 
84.2/101 260.4/164 


O(fix) 
37.5/50 


-20i| 
196.0/169 


71 +2i 
' i -13 

81.8/76 



Note — Terminologies and units arc same as table 2. 

The intensity of Fexxvi Ka is fixed at 30% of that of Fexxv Ka determined by the Suzaku fit. 
The other parameters are fixed to the best-fit Suzaku values. 
a Fine-tuning offset energy in units of eV. 



Table 2. The best-fit Suzaku (XIS) results of the Sgr B2 tematic errors, and refer to the Case 1 parameters in the 
region subsequent discussion. 

3.2.2. Dependence of the Line Flux on the Background 
Subtraction 

Parameter Terminology Value As we have already noted, the NXB of Chandra in 

the 2001 observation was unusually high. This unsta- 
ble (time-dependent) NXB is also often the problem in 
XMM-Newton analysis (Carter & Read 2007). The un- 
stable NXB and associated lines may affect the estimated 
line fluxes in the source spectra. To check these system- 
atic errors, we derived the count rates of MOS and PN 
for the source and background above the 10 keV band, 
where NXB is dominant. The count-rate ratios between 
the source and background were used as correction factors 
a for the background. The factors a for the XMM-Newton 
data in 2004 are relatively small, 1.0 and 1.4 for MOS and 
PN, respectively, while those for the data in 2001 are 2.1 
and 2.6 for MOS and PN, respectively. 

We multiplied the background spectra by the factors a, 
subtracted them from the source spectra, and then fitted 

them with the model expressed by equation (1) for the 

X 2 /dof 574.4/514"f pur cases (Case 1 to 4). The results are shown in table 

5. In this process, the CXB was also multiplied by the 
factor a. Therefore, we over(under)-subtracted the CXB 
by the factor a. The surface brightness of the CXB with 
a typical Galactic absorption of 6xl0 22 cm~ 2 is ~ 10~ 15 
ergs cm~ 2 s _1 arcmin -2 in the Sgr B2 region, which is 
only a few percent of the continuum flux in the region. 
Thus, possible over(under)-subtraction of the CXB can 
be neglected in the following discussion. 

From tables 4 and 5, we find that the fluxes of the 
6.40 keV and 6.67 keV lines with corrected background 
subtraction (table 5) are almost the same as those with 
normal background subtraction (table 4) within the sta- 
tistical errors, although the continuum flux is substan- 
tially different. Thus, in the following discussion based 
on the fluxes of the 6.40 keV and 6.67 keV lines, we can 
ignore possible uncertainty in the time-dependent back- 



— Continuum 



Absorption 




q+0.3 


Photon index 


r 


2 6 +0 - 2 

Z.D_q 2 


Power-law flux b 


± now 


c 9+0.4 


— Neutral iron lines — 




Center energy of Fe I-Ka c 


£^640 


6400+1 


Flux of Fel-Ka d 


^640 


1 14+0.04 


Center energy of Fel-K/3 C 


^706 


7058^ 


Flux of Fe I-K/3 d 


-P7O6 


ii +0 05 

u - li -0.07 


— Highly ionized 


iron lines — 




Center energy of Fe XXV-Ka c 


^667 


6665+^ 


Width of FeXXV-Ka e 


0667 




Flux of Fe XXV-Ka rf 


^667 




Center energy of Fe XXVI-Ka c 


^697 


6966113 


Flux of Fe XXVI-Ka d 


-^697 


n 09+0.05 

u - oz -0.04 



Note — Parentheses indicate the 90% confidence limit. 
a In unit of 10 23 H cm" 2 . 

b Observed flux of a power-law continuum in the 5—8 keV band in 
unit of 10 _ 4 photons cm~ 2 s _1 . 
c Line center energy in unit of eV. 

d Observed line flux in unit of 10 -4 photons cm -2 s — . 
e Line width in the Gaussian sigma in unit of eV. 
f The energy gap between Ka (6400 eV) and Kf3 (7058 eV) is 
fixed at the theoretical value (+658 eV) (Kaastra & Mewe 1993). 
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Table 4. Dependence of the best-fit flux on the continuum parameters for the Sgr B2 region 



MOS(2004) Fj 



PN(2004) 



MOS(2001) 



PN(2001) 



ACIS(2000) 



Parameter Casel Case2 
— XMM-Newton— 

+0.2 
^0.2 



Case3 



Case4 



pow 

Fq40 

^706 

^667 
P 

± pow 

-^640 

F7O6 

-F667 
P 

± pow 

^640 

F7O6 

-F667 
P 

± pow 
-^640 
^706 
-^667 



ACIS(2001) F l 



pow 

•^640 

^706 

-F667 
P 

± pow 
•^640 
F7O6 
■^667 



SIS(1994) F pow 

^640 
F706 

^ -^667 

Note — Same as table 3. 
Case 1 
Case 2 
Case 3 
Case 4 



4.8123 

1-0318:8 

"• lo -0.07 

qi+ - 04 

"• yi -0.08 

' - u -0.2 

5+0.05 
-0.07 

n 1 5+ - 05 

u - io -0.15 
-,+0.06 
-0.06 

a a+°- 4 

1 co+0.15 
i-O^-0.16 

+0.16 
-0.14 
+0.15 
0.16 
+0.3 
0.3 
+0.11 
0.12 
+0.09 
-0.09 
+0.11 
0.11 



1.06] 



0.89] 



0.16 

1.22 

9.4 

1.27 

0.17 

0.81 



4.9 

1 (u +0 - 08 

0.07 
0.07 
0.08 
0.04 
+0.2 
-0.2 
0.08 
0.05 
+0.05 
-0.06 
o+0. 05 
°-0.09 
+0.3 
0.3 
+0.12 
0.26 
+0.12 
-0.22 
+0.17 
0.19 
+0.3 
0.3 
+0.11 
0.11 
0.09 
-0.08 
0.10 
0.10 



0.10 
0.92 
6 
1.02 
0.06 

8.2 
1.53 
0.31 
1.16 
9.2 
1.19 
0.19 
0.72 



7.2 



-Chandra — 
6.6 



+0.5 



+0.. 



-0.5 

-1 n 7 +0.26 
1 -°'-0.26 

in +0 - 32 

u - lu -0.10 

,,+0.27 
-0.27 

6 7 +0 - 2 

"•'-0.2 
1+0.08 
-0.08 

17+ 009 

+0.08 
0.07 



1.14 H 



1.71" 



1.00 



■"-0.5 

1 %+ - 27 

l.OO_ 25 

12+ - 28 

U - 1Z -0.12 
q+0.27 
-0.27 
1+0.2 
-0.2 
h0.07 
-0.10 
q+0.12 
J -0.08 
F0.09 
-0.08 



1.16] 

6.0^ 

1.72] 

0.26^ 

1.02^ 



A S+ - 2 
4 - 8 -0.2 

1.07±g" 

+0. 
0. 
0. 
0. 
+0.2 
0.2 
0. 
0. 
+0. 
-0 

0. 
+0.4 
0.4 
0. 
0. 
0. 

1.10+°. 

9-6lg:S 
0. 

0. 

+0. 
-0 

0. 



0.12] 
0.95^ 
7.1] 
1.01] 

o.ioj 

0.85^ 
8.9^ 
1.43] 
0.16 H 



1.21] 
0.14] 
0.73^ 



7.1'-o.6 

1 oo+0. 27 
l.OO_ 24 

0.32 
-0.12 
0.14 
0.41 
+0.2 
-0.2 
+0.08 
0.08 
+0.08 
0.09 
0.08 

0.07 



0.12 

1.13 

6.5 

1.67 

0.21 

0.97 



0+0.2 

8 -0.2 
07+0-07 
u ' -0.07 
1 a+0.07 
iz -0.07 
qt -+0.07 
yo -0.07 
+0.2 
-0.2 
0.06 
0.06 
+0.04 
-0.05 
0.05 
0.05 
+0.3 
-0.3 
+0.14 
0.16 
+0.14 
0.14 
+0.14 
0.14 
+0.3 
0.3 
+0.11 
0.10 
+0.09 
0.08 
0.09 
0.10 



+0.4 

,,0+0.24 
48 -0.24 
9O+0.28 
zu -0.20 
on+0.24 
OU -0.24 
7+0.2 
' -0.2 
7fi +0.08 
'"-0.08 
+0.09 
-0.09 

n o+o.o7 

U8 -0.08 



31 





— ASCA 












5 


0+0.2 

°-0.2 


5 


c+0.2 

°-0.2 


5 


Q+0.2 

°-0.2 


5. 


1 


cq+0.18 
uo -0.17 


1 


"°-0.18 


1 


R9+0.14 
" z -0.19 


1. 





99+0.15 
zz -0.13 





on+0.17 
ou -0.15 





9z i+0.13 
z -0.15 


0. 


1 


90+0.18 
zu -0.18 


1 


14+O.I6 
1 -0.18 


1 


n7 +0.18 
"'-0.13 


1. 



65" 



1+0.2 

'-0.2 
+0.17 
0.10 
on+0.17 
°"-0.14 
1 c+0.15 
lo -0.16 



Nn an d r are free parameters (same as table 2 and 3) 
iVjj is a free parameter but T is fixted to the best-fit Suzaku value, 
r is a free parameter but TVji is fixted to the best-fit Suzaku value. 
TVh and T are fixted to the best-fit Suzaku value. 
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Table 5. Dependence of the XMM-Newton flux on the background-subtraction for the Sgr B2 region 



Parameter Cascl Case2 Casc3 Casc4 

XMM-Newton/PN(2004): Correction factor = 1.4 



p 

± pow 


c r+0.2 
°-°-0.2 


t; fi +0.2 
O.O_ 2 


c c+0.2 
O.O_ 2 


5.5_q 2 


-^640 


1 ri5 +006 

1 - uo -0.06 


1 n 7 +0.06 
1,u ' -0.06 




1.08 + '°-8| 

— U.UD 




06+ 05 

u - UD -0.05 


n 0fi +0 06 

u - uo -0.05 


07+ 05 


07+ 04 

u - u '-0.06 


^667 


qn +0 05 

u - au -0.06 


n OI +0.06 
u - yi -0.06 


Q2+ 06 

u - az -0.06 


Q2+ 06 

u - az -0.06 


X 2 /dof 


238.4/169 


240.9/170 


240.7/170 


241.1/171 


XMM- 


-Ncwton/MOS(2001): Correction factor 


= 2.1 


p 

± pow 


n+0.3 

°- y -0.3 


4 4+0-3 

^-0.3 


4 1+0-3 


4 6 +0 - 3 

^• u -0.3 


-^640 


1 40+O.I6 


1 ci +0.16 
i ' di -0.16 


1 co+0.16 
1 - oz -0.15 


1 46+ 015 


-P7O6 


21+ 015 


14+ - 14 


17+ 015 

' -0.14 


10+ 014 
u - lu -o.io 


^667 


1 90+O.I6 
1 - zo -0.16 


1 9 o+0.15 
l.zo_ 16 


1 oc+0.15 
J— '3-0.15 


1 17+0.15 

L - 1 1 -0.15 


X 2 /dof 


84.2/101 


92.6/102 


88.0/102 


94.7/103 



XMM-Newton/PN(2001): Correction factor = 2.6 

TP c n+O- 3 c 1+0.3 4 Q+0.3 r n+0.3 

r pow J - u -0.3 J - i -0.3 ^- y -0.3 d - u -0.3 

TP 1 98+0-12 1 oo+O.ll 1 or+0.03 1 o 4 +0.11 

i^706 0.24± ; 10 O^lto^g 0.24t°;o3 0.23l ; 09 

F 667 0.84±°;^ 0.87±g^ 0.92+™ 4 0.90±g : ~ 

X 2 /dof 260.4/164 267.5/165 267.1/165 267.9/166 

Note — Same as table 4. 



ground levels. 

3.3. The Constant 6.67 keV Line Flux 

The most serious problem when comparing the line 
fluxes between different satellites or different instruments 
is the reliability of the cross-calibration of their relative ef- 
ficiencies. The 6.67 keV line emission is due to the largely 
extended GCDX and a new SNR G 0.61+0.01, which con- 
tributes about 30% of the total 6.67 keV line emission in 
the Sgr B2 region. Therefore, the 6.67 keV line flux is 
time invariant and hence is a good cross-calibration line 
for the relative efficiencies of the Suzaku, XMM-Newton, 
Chandra, and ASCA observations. The constant flux hy- 
pothesis for the best-fit 6.67 keV line given in tables 2 
and 3 is rejected with % 2 /dof=56.3/8, which indicates 
that there are significant systematic errors in the detec- 
tion efficiencies of the ASCA, XMM-Newton, Chandra 
and Suzaku observations. 

The nominal systematic errors for these satellites are: 
13% for the ASCA SIS 3 , 5% for Chandra (Schwartz et al. 
2000), 10% for EPIC (also the MOS flux is 10-15% higher 
than the PN (Carter & Read 2007)) and 10% for Suzaku 
(Serlemitosos et al. 2007). 

The 6.67 keV line flux history obtained by adding these 
systematic errors to the statistical errors are plotted in 
figure 3. The % 2 /dof value with a constant flux model is 
then 4.9/8, which is acceptable at a 70% confidence level. 
Therefore, as expected, the 6.67 keV line flux is constant, 
if systematic errors are taken into account. 



> 
r- 

■•JD 

\6 




1995 



2000 
Time [year] 



2005 



Fig. 3. The time trend of 6.67 keV line fluxes in the Sgr 
B2 region. In this plot, each error bar indicates one sigma 
including the statistical and systematic uncertainties. The 
Suzaku data are indicated by triangles, XMM-Newton MOS 
data by circles, Chandra data by asterisks, and ASCA data 
by squares. The average flux is shown by the solid horizontal 
line. 

3.4- Time- Variable 6. 40 keV Line 

Unlike the 6.67 keV line, the 6.40 keV line emission is 
much more clumpy, hence it may not be time constant. To 
examine whether the 6.40 keV line from the Sgr B2 region 
is time variable or not, we performed a % 2 test using the 
same procedure as those of the 6.67 keV line. For the case 
3 See http://hcasarc.gsfc.nasa.gov/listserv/ascancws/msg00l43.html w ith no systematic errors, the constant flux hypothesis is 

rejected with x 2 /dof = 112.8/8. 

We next added the systematic errors in the same way as 
the 6.67 keV line case, and the constant flux hypothesis 
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Sgr B2 region [10~ 6 photons cm -2 s _l arcmins -2 ] 



1995 



2000 
Time [year] 



2005 



Fig. 4. The light curve of the corrected 6.40 keV line surface 
brightness in the Sgr B2 region. In this plot, each error bar 
indicates a 90% confidence limit. 

was still rejected with x 2 /dof = 17.9/8. Therefore, the 
6.40 keV line flux is time variable. 

Having confirmed that the 6.67 keV line is constant, we 
used it as a cross-calibrator. We introduced normaliza- 
tion factors for the detection efficiencies (c ± 8c) for the 
eight observations shown in table 6, where 5c is a statisti- 
cal error. These parameters arc the ratios of the best-fit 
6.67 keV line flux of each satellite to the averaged flux. 

The best-fit 6.40 keV line fluxes were divided by these 
factors, and are referred to as the corrected 6.40 keV line 
flux. Since the MOS and PN fluxes should be consistent, 
we averaged these fluxes and errors. Figure 4 shows the 
corrected 6.40 keV line fluxes in each observation period. 
The 6.40 keV line flux varies significantly with x 2 /dof = 
28.9/6. The flux of the 2005 observation is about 60% 
that of 2000. 

3. 5. Sub- structures in the Sgr B2 region 

In order to examine the morphology history of the 
6.40 keV line, we made a surface brightness map of iron 
lines for each observation. Since the energy resolutions of 
Chandra or ASCA were not sufficiently good to separate 
the 6.40 keV line from the 6.67 keV line, all the images 
are in the 6-7 keV band. For the continuum flux subtrac- 
tion from the 6-7 keV band flux, we assumed a power-law 
continuum of the photon index of 2.6, the best-fit Suzaku 
value, and estimated that the power-law continuum in the 
6-7 keV band is 65% of the 5-6 keV band flux. After 
subtracting the continuum emission, we made the surface 
brightness maps using a pixel unit of 50-arcscc square, 
similar to the PSF of ASCA, which had the worst PSF 
among the four satellites. 

Figures 5(a)-(d) show the surface brightness maps of 
the iron lines. The Suzaku and XMM- Newton images are 
relatively dim in the whole area compared to Chandra 
and ASCA, which is consistent with the spectral fitting 
result. In the Chandra image, a bright spot is found in 
the west, where a SNR candidate G 0.570-0.018 is lo- 
cated (Senda et al. 2002). We selected the brightest spot 



named M 0.66—0.02 (solid circle in figure 5) from its po- 
sition, which was termed the "Sgr B2 cloud" in previ- 
ous papers, and G 0.570—0.018 (dashed circle in figure 5). 
We then examined the 6.40 keV flux changes from these 
two spots. Since the Chandra data in 2001 did not have 
enough statistics for the analysis, we excluded this data 
hereafter. 

3.5.1. M 0.66-0.02 

We extracted the spectra of M0. 66-0. 02 (Sgr B2 cloud) 
from the solid circle with a radius of 3.2 arcmin in figure 
5 and subtracted the blank-sky data in the same detec- 
tor position. In the ASCA image, the peak position of 
M 0.66—0.02 is different from those in the others. It is not 
clear whether this shift is due to a relatively large error 
of the ASCA astrometry (no fine-tuning of the ASCA as- 
trometry was performed, see Section 2.4) or due to a real 
shift in the variable 6.40 keV line. Therefore, in the ASCA 
analysis, we selected two circles for M 0.66— 0.02, at the 
same as the position in the original coordinate of Suzaku, 
XMM-Newton, and Chandra (P) and at the center on the 
peak of the SIS image (Q). 

As described in Section 3.2, the dependences of the line 
fluxes on the continuum models and on the uncertainties 
of NXB (systematic errors) are not larger than the statis- 
tical errors of the Sgr B2 region. Since the statistical error 
of the line flux from M 0.66—0.02 is larger than that of the 
Sgr B2 region, these systematic errors can be neglected. 
We nevertheless attempted the same analysis as the Sgr 
B2 region and found no difference within different con- 
tinuum modeling and/or background subtractions. We 
therefore estimated the line fluxes using the continuum 
models of N H and T free. The ASCA/SIS fluxes of the 
6.40 keV line were 0.99^ for the spectrum region 'P', 
and 0.93^0 15 xl0~ 4 photons cm~ 2 s _1 for that of region 
'Q'. These fluxes are consistent within the statistical er- 
ror, hence we used the result of region 'P'. The best-fit 
spectra and parameters are shown in figure 6 and table 7, 
respectively. 

However, cross-calibration errors for the efficiencies be- 
tween the four satellites may not be neglected. We there- 
fore corrected the 6.40 keV fluxes using the correction fac- 
tors shown in table 6, which are derived from the 6.67 keV 
line flux in the Sgr B2 region (Section 3.5). The results 
are plotted in figure 7. This time history is similar to that 
of the Sgr B2 region (the whole region) (figure 4). The 
constant flux hypothesis is rejected with x 2 /dof=30.7/6. 
The flux in 2000 is 1.5 times higher than those in 2004 
and 2005, and declines in 2001. 

3.5.2. G 0.570-0.018 

G 0.570—0.018 is a diffuse source extending to about 10 
arcmin (Senda et al. 2002). We extracted the 6.40 keV 
flux from G 0.570—0.018 (the dashed circle with a radius 
of 2.5 arcmin in figure 5) in the same way as M 0.66—0.02. 
Since the spectra do not have enough statistics to derive 
Fel K/3 (7.06 keV) line, the intensity of K/3 was fixed to 
12.5% that of Ka (Kaastra & Mewe 1993). Using the 
same method as Section 3.5.1 except for the fixed line 
ratio, we fitted the Suzaku spectra with Nn and T free 
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Table 6. The normalization factors of the detection efficiency for the eight observations 



Suzaku XMM-Newton XMM-Newton Chandra Chandra ASCA 

XIS(2005) MOS(2004) PN(2004) MOS(2001) PN(2001) ACIS(2001) ACIS(2000) SIS(1994) 
1.05 ±0.05 0.91±0.08 0.89±0.06 1.25±0.14 0.83 ±0.11 1.12 ±0.29 0.99±0.07 1.19 ±0.18 

Note — Parentheses indicate in the 90% confidence limit. 

Normalization factors and their errors of the detection efficiency derived from the ratio of the best-fit 6.67 keV line flux in 
the Sgr B2 region obtained with each observation to the averaged flux. 
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Fig. 5. Surface brightness maps of iron lines obtained with (a) Suzaku XIS 2005, (b) XMM-Newton MOS and PN 2004, (c) Chandra 
ACIS-I 2000, and (d) ASCA SIS 1994. Pixel size is 50" X 50" in each case 
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Fig. 6. Spectra for M 0.66-0.02 obtained with Suzaku XIS FI sensors (a) and BI sensor (b) in 2005, XMM-Newton MOS (c) and 
PN (d) in 2004, (e) and (f) in 2001, Chandra ACIS-I (g) in 2000, and ASCA SIS (h) in 1994. The solid lines indicate the best-fit 
models. 



Table 7. The best-fit parameters of M 0.66-0.02 



Observatory 


Suzaku 


XMM-Newton 


XMM-Newton 


Chandra 


ASCA 


Detector 


XIS(2005) 


MOS (2004) 


PN(2004) 


MOS(2001) 


PN(2001) 


ACIS(2000) 


SIS(1994) 








— Continuum — 








r 

± pow 


r Q+0.5 
q+0.3 

2 o+ 01 

z - u -0.1 


°- -2.0 
q 9+O.5 
°- z -1.0 

1 6+ - 1 

i '°-0.1 


°- ' -1.1 

1 7+0.2 

i -'-0.3 

2 4+ 01 
z -^-o.i 

— Neutral 


2+ 24 

U.Z_ Q 2 

1+ 10 

q+0.2 

iron lines — 


9 r;+2.8 
Z '°-2.5 

1 2+ 11 
i - z -i.i 

2A±°i 


4 3+ 1 ' 1 
*-°-i.g 
1 9+0.5 

9 7+0-1 
z - ' -0.1 


9 q+0.2 
z '°-0.2 


^640 
-P7O6 

^667 


66+ 02 

06+ 04 

u - uo -0.04 

3fi+ 02 


2q+ 004 

u - za -0.04 


0.57t° Z 0.92+g;g 
0.02ig;g 0.051^ 
-highly ionized iron lines — 

n 91-+O.O3 n 4 9+0.09 
u - zo -0.03 u -^ z -0.09 


0-82t u:o ° 
11+ 008 

U-l-L-0 .07 

26+ 08 

u - ZD -0.07 


i-oo±85g 
0.0618:2 

25+ - 04 

u - zo -0.04 


n qq+013 

u ' yy -0.14 

o.i2l8:l? 

0.421°;!! 


^697 


u - lz -0.02 


0.10° 


0.08° 


0.14 a 


0.09 a 


0.08° 


0.14° 


X 2 /dof 


O(fix) 
348.4/354 


0(fix) 
36.9/40 


23l* 
71.3/68 


0(fix) 
35.3/38 


23±i 7 
16.6/27 


-19±§ 
142.1/144 


8ll 2 i 
36.1/37 



Note — Same as table 3. 

a The intensity of Fe XXVI Ka is fixed at 34% of that of Fe XXV Ka obtained with Suzaku. 
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(a) Suzaku XIS-FI 2005 



(b) Suzaku XIS-BI 2005 



(c) XMM-Newton MOS 2004 
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Fig. 8. Spectra for G 0.570—0.018 obtained with Suzaku XIS FI sensors (a) and BI sensor (b) in 2005, XMM-Newton MOS (c) and 
PN (d) in 2004, (e) and (f) in 2001, Chandra ACIS-I (g) in 2000, and ASCA SIS (h) in 1994. The solid lines indicate the best-fit 
models. 
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Fig. 7. The light curve of the 6.40 keV line surface brightness 
from M 0.66—0.02 (Sgr B2 cloud). In this plot, each error bar 
indicates a 90% confidence limit. 



model. We obtained 6415 eV as the line centroid of the 
6.40 keV line, which is 15 eV higher than the values ob- 
tained in the other region. We then refitted the spectra 
with the line centroids fixed at the values for the Sgr B2 
region (table 2) and found significant residuals in the 6.4- 
6.6 keV band. We added one Gaussian function in this 
energy band, and obtained an acceptable fit. The cen- 
troid of the additional line is 6505 eV, which is consistent 
with the value reported by Scnda et al. (2002). From the 
energy centroid, this line is likely to be due to FeXX-Ka. 
Assuming a thin-thermal plasma, the temperature of the 
most prominent emission of this line is 1.1 keV, a typical 
temperature of an SNR, which led Senda et al. (2002) to 
propose a new SNR in highly non-equilibrium ionization. 

Since the flux of the 6.51 keV line is due to highly ion- 
ized iron and the flux is likely to be time constant, the flux 
ratio among the highly ionized iron lines is fixed in the fol- 
lowing analysis. The best-fit spectra and parameters are 
shown in figure 8 and table 8, respectively 

We corrected the 6.40 keV fluxes using the correction 
factors as we did for M 0.66—0.02. Figure 9 shows the time 
trend of the 6.40 keV line from G 0.570-0.018, which in- 
dicates a clear variability in the 6.40 keV line emission. 
A constant flux hypothesis is marginally rejected with 
X 2 /dof = 13.2/6. 
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Table 8. The best-fit parameters of G 0.570-0.018 



Observatory 


Suzaku 


XMM-Newton 


XMM-Newton 


Chandra 


ASCA 




XTS( , 20n5 > ) 


MOS(2004) 


PN(2004) 


MOS(2001) 


PN(2001) 




kJTkJl t y JT 1 








— Continnum 








A H 


1.5±i-g 






1 3+ 1 - 4 






0-2121 


r 


2-6±giS 


4-i±i:S 




u -°-0.9 


2.o±i : « 


2 0+ 1 - 4 

z - u -0.6 


1 2 +i:s 
A - z -o.8 


± pow 


1 l +ai 
i - A -o.i 


i i +ai 

J - A -o.i 


1 3+ 01 
— Neutral 


1 q+o.i 

A ' y -0.1 

iron lines — 


1 2+ 01 


1 2+ 01 


q+ 01 


^640 


n ifi+ 03 


n 1 1 + 04 

u - AA -0.04 


n 14+003 

u - A ^-0.03 

— Ionized 


n n+0.06 

iron lines — 


n iq+006 

u - iy -0.06 


u - ZD -0.04 


12+ 005 


^667 


2Q+ 02 

u - ztl -0.02 


20+ 04 

u ' zu -0.04 


n 99+O.O3 
u ' zz -0.03 


n 1 r+0.06 


n 90+0. 06 
u - zo -0.06 


of;+o o3 

u - zo -0.03 


20+ 04 

u ' zu -0.02 


^651 


07+ 03 

U ' U '-0.02 


0.05 d 


0.06 d 


0.04 d 


0.05 d 


0.06 d 


0.05 d 



S Gain O(fix) O(fix) 35(fix) O(fix) 33(fix) -18 (fix) 71(fix) 

X 2 /dof 192.3/181 56.1/51 70.8/78 40.9/25 21.0/15 82.4/93 6.8/18 



Note — Same as table 3. 

d Fixed at the value and flux ratio against Fe xxv Kct line obtained with Suzaku. 



MO.66-0.02 [10 -6 photons cm -2 s _1 arcmins -2 ] 



1995 



2000 
Time [year] 



2005 



Fig. 9. The light curve of the 6.40 keV line surface brightness 
from G 0.570— 0.018. In this plot, each error bar indicates a 
90% confidence limit. 



Senda et al. (2002) mentioned that G 0.570-0.018 is 
a simple SNR candidate, but the 6.40 keV variabil- 
ity suggests that G 0.570—0.018 includes an XRN (see 
Discussion). In fact, Miyazaki & Tsuboi (2000) re- 
ported a flux peak of the CS (J=l-0) line emission at 
(I, b) = (0.°601, — 0.°025), which indicates the presence 
of a molecular cloud. The light curve of G 0.570—0.018 
is similar to that of M 0.66—0.02. Hence the distances of 
the two clouds from Sgr A* are the same. G 0.570—0.018 
may be located in front of M 0.66—0.02 with some off-sets 
of the line of sight. 

4. Discussion 

The 6.40 keV line is K-shell emission of neutral or low- 
ionization iron. Hence, since cool gases such as molecular 
clouds cannot emit X-ray photons, the 6.40 keV line is 
due to illumination by external sources of either X-rays 
or charged particles. The analyses of ASCA, Chandra 
and Suzaku (Koyama et al. 1996, Murakami et al. 2001b, 



Koyama et al. 2007b) revealed strong Fcl Ka,f3 lines and 
deep Fe K edge, which supports an X-ray origin. However, 
an electron origin cannot be excluded if the iron abun- 
dance is higher than the solar value by a factor of 3-4 
(Yusef-Zadeh et al. 2007) and the X-ray spot is at the 
back-side of the molecular clouds. 

The time variability of the 6.40 keV line in the Sgr B2 
region has already been reported by Koyama et al. (2008a) 
from the Chandra (2000) to the Suzaku (2005) observa- 
tions. A similar time variability has also been reported 
for the other cloud near Sgr A* by Muno et al. (2007) and 
Koyama et al. (2008b), which are the time variabilities of 
a few years span. In this paper, we have confirmed the 
time variability with a longer time scale (~10 years) from 
the Sgr B2 region. 

This time variability is also found in the sub-structures, 
M 0.66-0.02 and G 0.570-0.018. The sizes of these two 
sources are a few tens of light years, which is comparable 
to the variation time scale. These time variabilities are the 
clearest evidence for X-ray origin, because it is difficult to 
change the flux of any charged particle in a large region 
(a few 10 light-years) within a few 10 years. 

If the origin is X-ray irradiation, the observed 6.40 keV 
flux requires a bright source with a luminosity of more 
than 10 37 erg s _1 assuming the source is located near the 
Sgr B2 region in the distance of 8.5 kpc. Point sources in 
the GC region have been catalogued by many satellites, 
but no source has been brighter than 10 37 erg s _1 in a 
time scale of ten years. Thus, we conclude that the most 
probable source is the massive black hole Sgr A*. 

In this scenario, the light curves of the 6.40 keV line 
flux in the molecular clouds can be used as a tracer of 
the light curve of Sgr A* in the past, like a time-delayed 
echo. Since the light travel time between Sgr A* and the 
Sgr B2 region is ~300 years, we already proposed that Sgr 
A* was X-ray bright (~ 10 39 erg s _1 ) 300 years ago. Our 
new results indicate that the X-rays have been variable 
over a time scale of a few tens of years. The light curve 
of the 6.40 keV line is not the real history of Sgr A*. It 
is smeared out if the size of the reflection (fluorescence) 
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cloud is comparable or larger than the light-cross scale. 
Then, (1) even if the real light curve is pulse- like (delta- 
function), the observed light curve of the fluorescent line is 
smoothed on a scale of the light crossing time of the cloud, 
or (2) if the cloud is very small (like a point source), the 
observed light curve is the real history of the incident X- 
rays. The high-resolution Chandra image (Koyama et al. 
2008a) revealed that the fluorescent line emission has a 
complex structure with bright peak of a scale of less than 
1 arcmin (~6 light years). Therefore, a realistic situation 
lies between cases (1) and (2), so that the observed light 
curve is moderately smeared-out. This means that the 
observed variability amplitude is a lower limit of the actual 
one, although they do not differ largely from each other. 
Hereafter, we discuss the evolution of the Sgr A* outburst 
based on the observed light curve. The averaged decay 
rate of the 6.40 keV line is 60% in 10 years, hence the 
half-decay time is ~15 years. If this decay time has been 
constantly folded for 300 years, the X-ray luminosity after 
300 years becomes 10 33 ~ 34 erg s _1 , which is consistent 
with the observed value in the present. We hence propose 
that 300 years ago, Sgr A* experienced a giant flare and 
then entered a decay phase that continued to the present. 

5. Summary 

We summarize the results as follows: 

1. We have extensively studied the spectra of the 
Sgr B2 region, assuming various continuum mod- 
els, backgrounds and uncertainties of the relative 
efficiencies of Suzaku, XMM-Newton Chandra, and 
ASCA. We then found that the 6.40 kcV line fluxes 
are time variable. 

2. We found the time variability of the individual 
clouds, M 0.66-0.02 and G 0.570-0.018. 

3. The time histories of the 6.40 keV line flux are sim- 
ilar for all the sources, and exhibit their brightest 
peak in 2000, decline in 2001 and finally fall to 60% 
in 2005. 

4. The 6.40 kcV variability of the diffuse source 
G 0.570—0.018 suggests that this source is a can- 
didate for an XRN. 

5. The variability of the 6.40 keV line indicates that 
the GC black hole Sgr A* experienced a large flare 
300 years ago and entered a decay phase. 
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